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Motivated by anomalously large conductivity anisotropy in layered materials, we propose a simple 
model of randomly spaced potential barriers (mimicking stacking faults) with isotropic impurities 
in between the barriers. We solve this model both numerically and analytically, by utilizing an 
exact solution for the conductivity of a one-dimensional (ID) disordered system. In the absence of 
bulk disorder, electron motion in the out-of-plane direction is localized. Bulk disorder destroys ID 
localization. As a result, the out-of-plane conductivity is finite and scales linearly with the scattering 
rate by bulk impurities until planar and bulk disorder become comparable. The ac out-of-plane 
conductivity is of a manifestly non-Drude form, with a maximum at the frequency corresponding 
to the scattering rate by potential barriers. 

PACS numbers: 72.15.Rn, 73.20. Jc,73.21.-b 



It is usually the case that cleaner metals are better con- 
ductors. In the semiclassical, phase-incoherent regime of 
transport, this happens simply because stronger disorder 
means a shorter scattering time; in the phase-coherent 
regime, stronger disorder enhances Anderson localiza- 
tion thereby reducing the conductivity even further. It 
is also commonly believed that localization can be de- 
stroyed only by inelastic scattering. In this Letter, we 
propose and analyze a simple model with two types of 
disorder which defies these notions. We show that an 
increase in one type of disorder leads to a destruction 
of the Anderson-localized state and, consequently, to a 
increase in the conductivity in one direction. 

The model consists of planar barriers located at ran- 
dom spacings to each other and isotropic impurities dis- 
tributed randomly in between the barriers (see Fig. []}. 
This model is motivated by some well-known but hitherto 
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FIG. 1: (Color online) Left: a system of randomly spaced 
parallel potential barriers and randomly distributed isotropic 
impurities. Right: expected dependences of the in- and out- 
of-plane conductivities on bulk disorder. 

unexplained peculiarities of electron transport in layered 



conductors. In the band picture, the conductivity in a 
certain crystallographic direction scales with the inverse 
effective mass in this direction. In many cases, how- 
ever, the observed ratio of the in-plane and out-of-plane 
conductivities exceeds the (inverse) ratio of the effective 
masses by several orders of magnitude. A well-known 
case of such an anomaly is graphite, where the conduc- 
tivity ratio exceeds the mass ratio by 2-3 orders of magni- 
tude but other materials, e.g., NaCo2C>4 Q, cuprates 
[H, etc., also provide examples of this behavior. Stacking 
faults, e.g., "wrong" planes violating Bernal stacking of 
graphene sheets in graphite, have been proposed to be 
responsible for abnormally large conductivity anisotropy 
long time ago 0] ; however, little attention has been paid 
to localization of electrons by an array of faults. 

We consider a system of electrons with separable 
but otherwise arbitrary spectrum £r(fcii,fc z ) = + 
£z(kz), subject to two types of random potential: the 
ID potential of the barriers, U(z), and the 3D poten- 
tial of isotropic impurities, V (r). In the absence of 
bulk disorder, the in- and out-of-planes degrees of free- 
dom separate. Accordingly, the electron wave func- 
tion is factorized as V&mi,*:) = cp(r\\)x( z )> with 
x( z ) satisfying an effectively ID Schroedinger equation 

[s z (-id z ) + U(z)] X (z) = (E-e\\ where 

fcn is the (quasi) momentum along the planes. All states 
of such a system are localized in the z-direction by in- 
finitesimally weak disorder. Therefore, the dc conductiv- 
ity across the planes, a zz , is zero. On the other hand, 
since barriers do not affect the electron motion along the 
planes, the in-plane conductivity, mi, is infinite. Bulk 
disorder mixes the in- and out-of-planes degrees of free- 
dom, so that the separation of variables is no longer pos- 
sible. Therefore, ID localization in the z direction is 
destroyed, and a zz increases with bulk disorder, as long 
as it remains weaker than the planar one. When two 
disorders become comparable, a zz reaches a maximum 
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and decreases upon a further increase in bulk disorder in 
accord with the Drude formula [5J. At the same time, 
CT|| decreases monotonously with bulk disorder. A sketch 
of expected dependences of a zz and o\\ on 3D disorder 
is presented in Fig. [1] (right). In the rest of the paper, 
we confirm this simple picture both numerically, by cal- 
culating a zz in the Anderson model, and analytically, by 
exploiting the Berezinskii solution of the ID localization 
problem. 

Numerically, we study the Anderson model with 
nearest-neighbor hopping (set to unity to fix the energy 
scale) on a cubic lattice (of unit spacing) 
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Here, the on-site energy ei = tp\ + rji, and i = {i x ,iy,iz)- 
The first term, </>i, is the standard (bulk) disorder term 
which is chosen independently for each site in the inter- 
val (— W-b/2, Wb/2) with uniform probability. The sec- 
ond term, 7^ z , describing planar disorder, is chosen as 
—W with probability p and as W with probability 1 — p. 
For all results reported in this paper, p = 1/2. The sim- 
ulations are done at the energy equal to 0.1, to avoid 
the center of the band. We employ the recursive Green's 
function technique [|| with periodic boundary conditions 
in the directions transverse to the z-axis. The out-of- 
plane conductance G zz is equal to 2e 2 T/h, where T is 
the transmission coefficient between two wide leads. The 
simulations were performed for cubic samples of sizes L 
up to 35 lattice spacings. The bandwidths of planar dis- 
order W were chosen as 1,1.5,2,2.5 and 3, which cor- 
responds to localization lengths between roughly 2 and 
15 lattice spacings, in the absence of bulk disorder. The 
bandwidth of bulk disorder Wb ranged in between and 
18. We have averaged logG Z2 for 10 3 samples for each 
set of parameters. Crystalline anisotropy can be simply 
accounted for in simulations; however, the conductance 
is already anisotropic due to anisotropy of disorder even 
on a cubic lattice. 

Figure [5] shows G zz = exp{(lnG 2Z )} as a function of 
bulk disorder for several values of planar disorder. As 
expected, an increase in bulk disorder leads first to an in 
increase in G zz followed by a subsequent decrease. The 
position of the peak depends on planar disorder but is al- 
most independent of L. We checked that the conductance 
scales linearly with L for most of the range of parame- 
ters represented in Fig. 2, so that we are in the diffusive 
regime. Specifically, the diffusive regime begin when the 
conductance becomes larger than 2e 2 /h and continue up 
to the 3D Anderson transition (not shown in Fig. [J). 

Figure [3] shows the collapse of the data for the conduc- 
tivity, a zz — G zz /L, on a double- logarithmic plot. Three 
sets of curves corresponds to three values of planar dis- 
order: W — 1.5 (upper set), W = 2 (middle set) and 
W =2.5 (lower set). Within each set, the conductivity 
was computed for different values of L, as indicated in 




FIG. 2: (Color online) Out-of-plane conductance versus the 
bandwidth of bulk disorder Wb for a range of values of planar 
disorder W , as shown in the figure, and L = 30. 



the legend. The straight line has a slope equal to two. 
This scaling is confirmed by the analytic solution of the 
model, described below. 
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FIG. 3: (Color online) Out-of-plane conductivity versus the 
bandwidth of bulk disorder Wb on a double logarithmic scale 
for a range of system sizes, as shown in the figure, and three 
values of planar disorder: W = 1.5 (upper set), W = 2 (mid- 
dle set), and W = 2.5 (lower set). 



To solve the problem analytically, we adopt the delta- 
correlated forms for both types of disorder (U (z) U(0)) = 
-f z 5 (z) and (V (r) V (0)) = 7<S(r), and assume that 
bulk disorder is weaker than planar one, i.e., 1/r = 
2'nv^{Ep)^i <^ f/r 2 = 2-kv\(E, k\\)^ z , while planar dis- 
order is weak in a sense that EpT z ^ f . Here, is the 
3D density of states and v\ is the ID density of states at 
fixed value of k\ \ per one spin orientation. In the absence 
of bulk disorder, our problem reduces to the ID case 



with the velocity v z = \de(k\\,k z )/dk z 



where 



=fe*F(fc||)' 

k Z F(k\\) is a positive root of the equation e{k\\,k z ) = Ep, 
and the scattering time t z being functions of fen. The re- 
sult for the ac conductivity of a strictly ID disordered 
system, surmised first by Mott 0] and derived rigorously 
by Berezinskii Q reads 



16e 2 v z r z 



i((3)ut z + 2t 2 lo 2 \ti 2 (wt,)] ,(2) 



for ljt z <C 1. (The numerical coefficient in the imaginary 
part was corrected in Refs. [UHH). The out-of-plane con- 
ductivity of a 3D sample with V — is obtained from Eq. 
(T2]) by summing over fen: a zz (u>) = J d 2 k\\o 1D (a>) /(2tt) 2 . 
As expected, a zz (0) = 0. 

In the presence of both types of disorder, a zz is given 
by the Kubo formula 



ffa(w) = 0S I dZ ' 



fcn,fc' 



x ((«,^(fc||, z; fcf,, z'KSf (fcf, , z'; fe,,, z)) p ) b , (3) 

where = g R( - A \k\\, z;^, z'; Ep ± w/2) is an ex- 

act retarded (advanced) electron Green's function in the 
mixed kn—z representation for a given disorder realiza- 
tion, A is the sample area in the lateral direction, and 
(. . . }b, P denotes averaging over bulk and planar disor- 
ders, correspondingly. The diagram for a zz is shown in 
Fig. Q]on the left. To leading order in 7, the conductivity 
VzJ averaged over bulk disorder is given by the sum of 
the two diagrams in the first row of Fig. |4] where thick 
solid lines denote Green's functions in the absence of bulk 
disorder, Gp^ (z, z'\ fen; E), and zigzags denote the cor- 
relation function of bulk disorder. There are no vertex 
corrections for the case of delta-correlated bulk disor- 
der. The first (second) diagram in the first row of Fig. [4] 
is obtained by replacing the exact Green's function by 
lf tu p u Gp A \z,Zi)k\\)Gv { A) {z x ,zi;p\\)Qp { A '{zuz';k\\). 

Subsequent averaging over planar disorder is simpli- 
fied dramatically by noticing that the effective energies 
E — £||(fc||) of the Green's functions depend on a 
particular value of fen. For short-range bulk disorder, 
the momentum pji of the Green's function below the 
zigzag line differs considerably from the momentum fen 
in the rest of the diagram. This means that the typical 
difference of corresponding energies is of order Ep, 
i.e., much greater than 1/t 2 . In this situation, one can 
safely neglect correlations between the Green's functions 
with different momenta and average Gp(z\, Z\;p\\, E) 
over planar disorder independently from the rest of 
the diagram. As a result, we arrive at the diagrams 
in the second row of Fig. [4] where thick dashed lines 
denote the Green's function averaged over planar 
disorder. For weak planar disorder (EpT z 3> 1), 




FIG. 4: (Color online) Diagrams for the out-of-plane conduc- 
tivity to leading order in bulk disorder. Thin lines: exact 
Green's functions in the presence of both types of disorder; 
thick solid lines: Green's functions in the presence of planar 
disorder only; thick dashed lines: Green's functions averaged 
over planar disorder; zigzag: correlator of bulk disorder; solid 
and dashed brackets: averaging over planar and bulk disorder, 
correspondingly. 



this Green's function 
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e,(A*)-e|| (fe||) it/ 2 **] /2tt 
the corresponding self-energy insertion reduces 
to E R(A) (z,z';k {l ;E) = ^(i/2r)6{z - z'). Ex- 
panding G R ' A over the basis exact of eigenstates 
of the ID problem, we reduce the convolution 
of two Green's functions, sharing the point z\, 

to L^p' A {^^A\^ E )^' A { z ^ z '-A^ E ) 

~~§e'3p ,A ( Z ' 2: ';£||;- E )- Consequently, 0$ (u) is 
obtained from the exact ID result via 



r (i) 



, / d 2 fc|, da 1D (w) 



(27T) 



duj 



(4) 



To obtain the dc conductivity, one needs to differentiate 
only the imaginary part of Eq. (J2J) . This gives 



a£(0) = 2e 2 v 3 (Ep)D z 



where 



16C(3) 



\VUk\\))\\ ^,max 



(5) 



(6) 



z,max and ti 2jlnax denote the maximum values of 
= v z t z and u z (fei|) attained for fen =0 and 



Here. I 

lz{k\\) 

(/(fell)),, = (A^{Ep)) 1 /d 2 fe||^i(^,fe||)/(fc||). The 
diffusion coefficient D zz is proportional to the ratio of the 
square of the localization length in the ID system to the 
bulk scattering time. Numerically, we have found that 
a zz (0) scales as the square of the bulk disorder band- 
width. This is confirmed by our analytic result since 
Ozz (0) oc 1/t oc 7 cx B 2 in the Born approximation. 

Equation (J6j) allows for a simple physical interpreta- 
tion. Bulk scattering weakly couples ID channels of 
localized electrons with different fen. Each scattering 
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event results in a random displacement of order l z in the 
z-direction and results in diffusion with the coefficient 
D zz ~ Iz/t. Notice that bulk disorder acts very similarly 
to the electron-phonon (e-ph) interaction in a strictly ID 
system, where a 1D (0) oc l/r _ p h [9]. The difference be- 
tween the two cases is that a 1D (0) scales with l/r _ p h 
only at temperatures higher than the single-level spacing 
within the localization length, i.e, for Tt z 3> 1, while at 
lower temperatures a 1D is of the hopping form. The con- 
dition Tt z ^> 1 allows one to neglect correlations between 
the Green's functions in the self-energy insertions and in 
the rest of the diagram. In our case, these correlations 
can be always neglected for short-range bulk disorder, 
i.e., in contrast to phonon-activated transport, there is 
no "hopping" regime for disorder-activated transport. 

Coming back to the issue of anomalously large con- 
ductivity anisotropy, it is easy to show that the in- 
plane conductivity is given by the usual Drude formula 
<j Q( 9 = 5 a pG 2 vj{EF)(v a Vj3)] i ] i T . Then the conductivity ra- 
tio can be estimated as 

cr\\/a zz ~ ((«ij)||/<^>||) (r/T Zimax ) 2 . (7) 

As an example, we consider the case of graphite with 



&\\/&zz — 10 4 at low temperatures. A realistic band 
structure model of graphite [l[ gives {vf<} / {v 2 ) ~ 140, 



11] 



thus t z /t ~ 0.12. Taking r = 0.5 x 10~ 12 s from Ref. 
and estimating (v 2 ) 1 ^ 2 ~2x 10 6 cm/s, we obtain for the 
mean free path due to planar disorder (stalking faults) 
m l z ~ 120 A. This means that stalking faults are sepa- 
rated by about a hundred perfect planes, which is quite 
a realistic assumption. 

Summing up higher-order diagrams with self-energy in- 
sertions due to bulk disorder amounts to replacing the ex- 
act Green's functions in Eq. © by Q R{ - A \z : z'\ % E ± 
w/2±^), which can be viewed as functions of a complex 
frequency. One can verify that all intermediate steps in 
Refs. [H and [1(3] are valid for complex lu as well. There- 
fore, the general result for the conductivity of our model 
is obtained from the Berezinskii's solution as 



<? zz (lu) 



(2^V 



-l-D 



LU + - 



(8) 



To lowest order in 1/r, Eq. §8§ reduces back to Eq. (@|. 
Within the logarithmic accuracy of the original Berezin- 
skii's formula, we obtain 



Recr zz (w) 
lma zz (u>) 



= 2eV 3 



16/2 



C(3) 



^ (lu 2 t 2 - l) In 2 \ lu 2 t 2 + Tr 



-2e v^lut ( 



/ml 



C(3)-^ln 2 L 2 r 2 + J| 



(9a) 
(9b) 



These formulas are valid for an arbitrary value out but 
only for lut z <C 1 and t z /t <C 1. From Eq. (|5a|) . we 
see that Kea zz (uj) is almost constant for to <C <-u cr = 

1/2 

1/ \Tz,raaxT) and increases with w in a Mott way, as 
w 2 ln 2 ti;, for lo ^> co cr . At higher frequencies, uj l/r z , 
o zz (ix>) can be found perturbativcly in l/r 2 : the lead- 
ing order result is simply a Drude formula a zz {oj) oc 
1/lu 2 t z —i/oj. Therefore, both Kea zz (u>) and — lma zz (u>) 
have maxima at u> ~ 1/t z . Thus, although bulk disor- 
der destroys localization at lu — 0, the resulting state still 
has properties interpolating between those of a metal and 
an Anderson insulator. This prediction is amenable to a 
direct experimental verification. 

Finally, we notice that the predictions of our model are 
equally well applicable to a two-dimensional (2D) case, 
e.g, for line barriers crossing the plane. Such a system 
can be realized in a 2D electron gas with an array of 
randomly spaced stripe-like gates. 

In conclusion, we have shown that a system with two 
types of disorder-randomly spaced planar barriers and 
bulk impurities- exhibits quite unusual transport prop- 



erties. In the absence of bulk disorder, it behaves as a 
ID insulator in the out-of-plane direction and as an ideal 
metal in the in-plane direction. Bulk disorder renders 
both conductivities finite; however, a zz increases with 
bulk disorder until two disorders become comparable. 
For weak bulk disorder, the ratio of the conductivities 
may exceed the ratio of the effective masses by orders of 
magnitude. The ac out-of-plane conductivity has a mani- 
festly non-Drude frequency dependence with a maximum 
at intermediate frequencies. 
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